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ABSTRACT 

A new technique is presented for thermographic detection of flaws in materials and structures by performing temperature 
measurements with fiber Bragg gratings. Individual optical fibers with multiple Bragg gratings employed as surface 
temperature sensors were bonded to the surfaces of structures with subsurface defects or thickness variations. Both 
during and following the application of a thermal heat flux to the surface, the individual Bragg grating sensors measured 
the temporal and spatial temperature variations. The investigated structures included a 10-ply composite specimen with 
subsurface delaminations of various sizes and depths. The data obtained from grating sensors were further analyzed with 
thermal modeling to reveal particular characteristics of the interested areas. These results were found to be consistent 
with those from conventional thermography techniques. Limitations of the technique were investigated using both 
experimental and numerical simulation techniques. Methods for performing in-situ structural health monitoring are 
discussed. 

Keywords: fiber Bragg grating, temperature sensor, thermal health monitoring, thermography, composite, delamination 


1. INTRODUCTION 


Fiber optic sensors have been extensively studied and proposed for temperature and strain sensing in heath monitoring 
systems of aerospace structures and materials as well as many other applications. 1 " 3 Compared to other sensors, Fiber 
optic sensors have the advantages of being light-weight and flexible, and requiring simpler wiring especially for 
distributed sensing. For an extensive heath monitoring system distributed (multiple) sensors are an important 
requirement. Most efforts have focused on using distributed fiber optic strain sensing systems, with limited consideration 
given to temperature sensing systems. For fiber optic temperature sensing, the conventional techniques using Raman 
scattering can detect the average temperature over a long distance within a single fiber. However the spatial resolutions 
for temperature readings-averaging over a length of fiber-are typically about one meter. Therefore they are not suitable 
for measuring the thermal responses at specific locations of the investigated materials. 

Another more attractive alternative is using fiber Bragg gratings (FBGs) for distributed temperature sensing (DTS). 
FBGs written in an optical fiber have a physical length of only a few millimeters. This size is nearly a point sensor. 
However, if multiple FBGs are written at different (multiplexing) wavelengths and read with a time domain 
demodulation system, the total number of these FBGs in a singe fiber is limited due to the finite bandwidth of the laser in 
the detector system. To solve the above problems of spatial resolution and limited sensors, recently a new technique for 
fiber optic DTS has been developed using low reflectivity FBGs. Under this technique low reflectivity FBGs were 
written at the same nominal wavelength and read with a frequency domain demodulation system. 4,5 This technique 
allows a single fiber to contain hundreds of FBGs, employed as temperature and strain sensors. 

An effective fiber optic DTS system has some advantages compared to a distributed strain sensing system for structure 
health monitoring. For a strain sensing system, in general, it is required to mechanically excite (or stress) the investigated 
structure or material. 6,7 For some cases this is impractical, if not impossible. Obviously using a DTS system requires the 
investigated material to be thermally excited. However the temperature change of the material from the required thermal 
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excitation can be only a few degrees or no more than tens of degrees. Such temperature variation is comparable with the 
temperature changes excited for conventional thermography techniques. 

Conventional thermographic techniques in general utilize a flash or quartz lamp as a heating source and an infrared 
imager to detect the thermal response of the investigated material. The heating pulse duration typically ranges from a 
fraction of a second to a few seconds. The IR imager contains an array of several hundred by hundred detectors. These 
techniques are capable of large area inspection of aerospace structures and materials for their reliability and safety. One 
of the particularly important areas is the inspection of graphite fiber reinforced composite materials because they are 
being increasingly used as primary structures due to their high stiffness and strength to weight ratio. Of particular 
interest is the detection of delaminations that can appreciably reduce the compressive strength of a composite. Recently 
Winfree et al. have successfully developed a physical model of two-layered systems for accurate reduction of the 
temporal thermal response of a composite with fabricated delaminations to the depths of the delaminations. 8,9 

In this paper we propose a new technique using a single optical fiber with multiple FBGs for the thermographic 
detection of flaws in materials and structures. The investigated structures included a 10-ply composite specimen with 
fabricated delaminations of various sizes and depths, similar to the one mentioned above. The optical fiber was bonded 
to the surface of the investigated composite. Both during and following the application of a thermal heat flux to the 
surface, the individual Bragg grating sensors measured the temporal and spatial temperature variations. The data 
obtained from individual FBGs were analyzed with thermal modeling to reveal particular characteristics of the interested 
areas. 


2. THEORY 


2.1 FBG as a temperature sensor 

In general a fiber Bragg grating can be characterized by its Bragg wavelength, which is the center wavelength of the 
light reflected from the grating. The Bragg wavelength is given as 


A. B -2n e ffA, (1) 

where n e ff is the effective refractive index of the fiber core and A the grating period. For a fiber Bragg grating bonded 
onto or embedded in a polymeric substrate, a change in the temperature causes a change in the grating period due to not 
only the thermal expansion of the fiber but also the strain induced by thermal expansion of the substrate. In addition, the 
refractive index of the fiber core changes because of the thermo-optic effect. Combining all the above effects, the shift 
in the Bragg wavelength due to a temperature change, AT, is given as 

bX B Ab = ( 1 - Pe)8 / // + 6n e ff/n e f f , (2) 

where p c is the photo elastic constant of the optical fiber and SI /I is the thermally induced strain of the fiber. In general 
the photoelastic constant, the fractional index change, and the thermal expansion coefficients of the fiber and the 
substrate are temperature dependent and nonlinear, especially at low temperatures. Flowever, for a finite temperature 
change, especially for the temperatures around and above the room temperature, Eq. (2) can be rewritten as a linear form 


6X B A B = [(l-pe)a p +i]AT, (3) 

where § = (l/n e ff)( <?n e ff I dT), is the thermo-optic coefficient of the fiber and a p is the thermal expansion coefficient of 
the substrate. Eq. (3) takes into account the conditions when the thermal expansion coefficient and the physical 
dimension of the substrate are much greater than those of the optical fiber. This simple equation allows the FBG to 
perform as a temperature sensor of the substrate. 



2.2 Determination of delamination depth from the theoretical thermal response of two-layered system 

Two steps are required to determine the delamination depth from the thermal response of a composite. The first step is 
the development of a model for the thermal diffusion in the specimen. For this case, a one-dimensional model is 
proposed in the following section. After development of a model, the thermal data must be compared to different model 
responses to determine which model parameters give the best approximation of the data. 

Real delaminations in composites are air gaps between two layers in the composite. It is well known that for long times, 
the response of the temperature over the delamination is not characteristic of a single layer response. The typical 
explanation for this effect is that lateral heat flow around the delamination diffuses the heat to the backside of the 
delamination. However if the thickness of the air gap is such that the heat flow across the gap is much larger than the 
heat flow around the delamination, the gap’s thermal resistance dominates the time response of the delamination. Note 
that if the thermal resistance of the delamination is sufficiently small, the time response at the front surface may not be 
significantly different than a single layer response. 

A simple analytic solution does not exist for the one-dimensional heat flow in a multilayered material. A solution does 
however exist in Laplace space for two layers of thickness I\ and l 2 coupled by an intermediate contact resistance (R). 
Since the configuration of interest is a composite with a delamination, the first and second layers are assumed to have the 
same thermal conductivity (K) and diffusivity (k). For the surface with the incident heating, the Laplace transform of the 
temperature response is given by 9 

T = f(p) (sinh(7 1 q) sinh(/ 2 q) + cosh (/ x q) (cosh(/ 2 q) + Kg R sinh(/ 2 q))) 

1 Kq (sinh(/, q) cosh(/ 2 q)+ (cosh(/j q) + KqR sinh(/ 1 q)) sinh(/ 2 q)) 


where q is \ p/K ,p is the coordinate in Laplace space and 

f(p) = foV-e- pt °)/p 

where fo is the energy per area from the heating source and t 0 is the heating duration. 

The parameter of interest in this study is the depth of a delamination. The depth of a delamination as measured from the 
front surface is 1\ in Eq. (4). Since the depth always appears with diffusivity in the specimen, it is not possible to obtain 
a value of the depth from the thermal response without an independent measurement of diffusivity. It is however 
possible to determine the relative depth of a delamination, if there is a clear region of the composite where no 
delamination exists. Representing the relative depth of the delamination as r, (/i = r l and 1 2 = (1 - r) I), Eq. (4) can be 
rewritten as 


T f (p) 


a(p)( 2cosh(^) + J^R K (sinh(J^) - sinh(^-2 

(cosh - cosh( ^ - 2 t 



(5) 


where a(p) is j[p)p IK, R K is RK //, and K/ is the diffusivity of the specimen divided by the full thickness squared. K/ is 
estimated from a fit of the region of the specimen where there is no delamination. 

With this parameterization of the equation, it is possible to find the relative depth of a delamination. Before finding the 
delamination depth, the characteristic K/ is determined from either a single-side or two-sided diffusivity measurement for 
the whole specimen. 9 This enables both an estimate for k / and a good image for identification of delaminations. The 
characteristic K/ is estimated by taking the median of the values for k /. 

For a longer duration heating from a quartz lamp, the thermal responses at sensor locations can be calculated by 
convolving these impulse responses. 



3. EXPERIMENT 


The low reflectivity FBGs used in this research were written in situ, into the optical fiber drawn using the NASA 
Langley optical fiber draw tower. A pulsed KrF-excimer laser of 248 nm and a Talbot interferometer arrangement 
mounted on the tower were used to write the FBGs. The interferometer consisted of a phase mask functioning as a beam 
splitter and a pair of mirrors used to recombine the split beams to form an interference pattern. An aperture placed in the 
laser beam path was used to control the grating length and spatial profile." FBG lengths for the present study were 
nominally 5 mm. The grating pitch written into the fiber could be adjusted by changing the relative angle of the two 
mirrors. These single-mode fibers were drawn from commercially available germanium-doped preforms of high 
numerical apertures. The drawn fibers with FBGs were coated with polyimide to thicknesses ranging from 11 to 16 
micrometers, then ink-marked to show grating locations. 

These low reflectivity FBGs (single grating reflectivity smaller than a few tenths of one percent) were interrogated using 
a frequency domain demodulation system shown in Fig. 1. In this system, the fiber coupler Ci and a pair of Faraday 
rotation mirrors (FRMs) form an in-fiber interferometer with an optical path difference of 2n e ffL 0 , where n e ff is the 
effective refractive index of the fiber core and L 0 the length of the reference cavity. The signals are driven by the tuning 
of the laser and detected at the photo-detector D;. They are used to trigger the sampling of signal at D 2 , which is the 
output of another in-fiber interferometer formed with the fiber coupler C 2 , a broadband reflector, and a particular fiber 
Bragg grating at a distance of Lj. If there is a series of low reflectivity Bragg gratings written at the same wavelength on 
a single fiber at different locations, the reflected signals from each grating are superimposed and detected at D 2 . The 
detected signals are further processed (fast-Fourier- transformed) to obtain the spatial spectrum of all gratings, which 
displays the physical profiles of the gratings at different locations. The spatial spectrum of a particular grating can then 
be windowed and inverse-fast-Fourier-transformed to get its own wavelength spectrum. Fig. 2(a) shows the spatial 
spectrum of two FBGs in a fiber with multiple gratings. The gratings have a physical length of about 5 mm and a 
separation distance of 10 cm. Two lines in Fig. 2(a) indicate the window of the inverse FFT for the second FBG. The 
calculated wavelength spectrum of the second FBG gives a center wavelength of 1551.3 nm, as shown in Fig. 2(b). For 
the distributed thermal sensing, the center wavelength shifts of each FBG were measured and calculated using the 
process described above. 



broadband FBGS 
reflector 


Fig. 1. Schematic diagram of a frequency domain demodulation system. Items Co, Ci, and C 2 are fiber couplers. “X” 
indicates that the unused port is terminated. Items Dj and D 2 are detectors and FRMs are Faraday rotation mirrors. 


In this study, an optical fiber with 32 FBGs was bonded to an investigated composite specimen with Capton® tape, as 
shown in Fig 3 (a). This bonding technique allowed the optical fiber to be bonded and taken off the sample quickly and 
neatly. At the same time it also allowed the fiber to sense the thermal strain for a temperature change up to 50° C 
without imminent slipping. 
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Fig. 2. (a) The spatial spectrum of two FBGs with a physical length of about 5 mm and two lines indicating the window of 
the inverse FFT for the second FBG. (b) The calculated wavelength spectrum of the second FBG with a center 
wavelength of 1551.3 nm. 


The specimen used for testing the viability of this measurement technique was a composite specimen with 20 
delaminations at specified depths. The 10 ply quasi-isotropic composite panel with a lay-up of [0,45, 90, -45,0/, 45, 90,- 
45/, 0,90/] was 31.75 x 31.75 centimeters and 0.19 centimeters thick (Fig. 3(a)). The delamination defect areas were 
squares with sizes of 3. 8x3. 8, 2. 5x2. 5, 1.9x1. 9, and 1.3x1. 3 square centimeters. The defects were buried at depths of 10, 
20, 30, 40, and 50 percent of the total thickness. A schematic of the defect layout and FBG sensor numbers is shown in 
Fig. 3 (b). The optical fiber bonded to the specimen had FBGs with a separation (center-to-center) distance of 5.08 
centimeters, which was the same as the separation distance of the squares in the same row. One FBG was placed at the 
center of each square. Five FBGs numbered 14-19 were outside the squares for comparison. 



Fig. 3. (a) A single optical fiber with multiple Bragg grating sensors bonded onto the surface of a 10-ply-composite, (b) The 
fabricated delaminations of various sizes arranged in columns between plies 1-2 (A), 2-3 (B), 3-4 (C), 4-5 (D), 5-6 (E). 


A quartz lamp was used to heat the front surface of the specimen. The heating durations were typically about 10-15 
seconds. The grating data acquisition was performed during and following the application of the heat flux to the surface 
at a rapidity rate of about 2.85 Hz. An IR camera was also used to record the temperatures of the surface separately. 
The experimental setup is shown in Fig. 4. 







Fig. 4. The experimental setup of the fiber optic thermal health monitoring system with a quartz lamp as the heating source 
and a conventional 1R camera for comparison. 


4. RESULTS AND DISCUSSION 

The time evolution of the temperature of the composite during the heating and cooling is a function of its three- 
dimensional geometry and thermal properties. For a thin composite the initial temperature response is initially a result of 
heat diffusion normal to the surface and is approximately one-dimensional and given by Eq. (5). As time progresses, 
first the in-plane heat diffusion, then the convective heat loss at the surfaces dominates the temperature response. The 
expected temperature response for the composite with delaminations during and after a quartz lamp heating pulse was 
estimated using a finite element routine. 

Common literature values of thermal properties of the composite were used. The thermal conductivity perpendicular to 
the fibers was assumed to be 0.97 W m^K’ 1 and 4.85 W nf'K' 1 along the fiber. The specific heat and density were 
assumed to be 1270 Jlsg^K" 1 and 1490 kgm" 1 respectively. A heat transfer coefficient of 8.063 W m" 2 K _1 for convection 
loss at the front and back surface was used, based on the exponential decay in temperature at late times as measured by a 
fiber optic sensor at the center of the panel. For delaminations, 0.025 cm air gaps were placed at depths that 
corresponded to the interface between plies at positions corresponding to the squares in the schematic of the composite 
shown in Fig. 3(b). 

In the simulation, first a short heating pulse of 10 ms duration was applied to the outer surface of the composite. The 
temperatures at the front and back surface were calculated at 10 ms intervals for the first 3.5 seconds, then at 3.5 second 
intervals for a total of 193 seconds. The response for quartz lamp heating was calculated at sensor locations by 
convolving the response from this “impulse” response, with a pulse shape that is estimated from the fiber optic sensor in 
the center of the panel to be a 2 second linear ramp up to a constant flux, with 13.5 second hold and a subsequent 2 
second linear ramp down. Using this pulse shape, a model response at a sensor location without delamination was 
calculated as shown in Fig. 5. 

As described in Section 3, the thermal responses of each FBG in the fiber were calculated from the measured center 
wavelength shifts of the FBG in a temporal sequence. Fig. 5 also shows the temporal thermal response of the FBG 
Sensor 16, located in an area without delamination. The relatively good agreement can be seen between the fiber optic 
sensor and the model response. 
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Fig. 5. Comparison of finite element simulation to fiber optic measured temperature response at center of composite in 
region with no delamination. 
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Fig. 6. Comparison of finite element simulation to fiber optic measured temperature response at center of composite in 
region with a delamination. 

Similarly, shown in Fig. 6 are the simulation and the measured fiber optic temperature response for an area with a 
delamination depth of 10 percent of the total thickness. Both cases show good agreement between the simulation and 


the experimental data during and immediately following the heat pulse. At later times the temperatures measured by the 
fiber optic sensor significantly deviates from the simulation results. This deviation is most likely the result of uneven 
heating, which is not incorporated into the simulation. 

By neglecting the thermal response of the later cooling stage and comparing the measured fiber optic data with the 
model, the depths of delaminations can be estimated. Fig. 7 shows the estimated depths of the delaminations from the 
fiber optic measurements for all of the 3.8 cm x 3.8 cm squares (no. 1, 3, 5, 29, and 31, as shown in Fig. 3 (b)). The 
solid line represents the model values of the depths. The first one (10 % of the thickness) has the highest degree of 
agreement (0.91, the ratio of estimated depth/ specification depth). The rest have lesser degrees of agreement, ranging 
from 0.55 to 0.72. These results are consistent with those from the conventional thermography technique applied to the 
same sample, 9 but with greater errors, especially for the deeper delaminations. 

From the current set of measurements, the depth of the smaller sizes of delaminations is more difficult to determine. 
Based on the simulations, the physical length of FBGs ( typically 5 mm) does not seem to be a limiting factor for these 
measurements. The depth measurement for the smaller delaminations should improve if the heat pulse duration is 
reduced and the frequency of the temperature measurement is increased. These changes will reduce the effects of 
uneven heating and result in the thermal response more closely approximated by the one dimensional model . 



Fig. 7. Comparison of the estimated depths of the delaminations measured by the optic sensors to the model values for all of 
the 3.8 cm x 3.8 cm squares. 


5. CONCLUSION 


We have developed a new technique of using distibuted FBG sensors for thermographic detection of flaws in materials 
and structures. Individual fibers with multiple FBGs employed as temperature sensors can be bonded to the surfaces of 
structures or embedded in the structures. By applying a thermal heat flux to the surface of the investigated structures, 
the individual Bragg grating sensors successfully measured the temporal and spatial temperature variations on the 
surface. These thermal responses were consistent with the results from conventional thermographic techniques. 9 We 
have also applied this technique on a 10-ply composite with fabricated delaminations of various sizes and depths. With 
different degrees of success, it has detected the square delaminations with widths of 3.8 cm at depths from 10 % to 50 % 


of the total thickness. Future efforts will focus on developing the technique with a faster detection system and assessing 
its potential for performing thermal health monitoring of aerospace structures and materials. 
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